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ABSTRACT 


A mathematical model of two Center Pivot sprinkler irrigation sys- 
tems (Circle Master and Valley 1160), operating under local climatic con- 
ditions, was developed to examine the effect of wind on the uniformity of 
water distribution and to evaluate the systems performance under actu- 
al windy conditions. 

The model consists of two parts. In the first part, the sets of 
equations of motion were solved for known initial conditions by means 
of a Continuous Simulation Modelling Program (CSMP). The simulation 
model described the trajectories of the particles in flight which per- 
mitted defining the geometry of the basic wetted patterms when the 
water particles struck the ground surface. The second part of the mod- 
el simulated the relative performance of two operating C.P. systems in 
relation to the rows of measuring receptacles used in field measure- 
ments. This involved determining the actual wetted patterns from the 
basic wetted patterns, Simulation of distribution patterns inside the 
actual wetted patterns and simulation of climatic factors and systems 
characteristics. The Fortran IV computer language was used to program 
the necessary formulation simulating the two operating C.P. systems. 

Performance of the simulated model was evaluated by a direct com 
parison of field and theoretical simulated application depths and by 
comparison of the theoretical simulated coefficients of uniformity with 
those measured. On the whole, theoretical depths deviated approximate- 
ly 14 percent from the field measurements. The spread of average 
deviations varied from 1l to 18 percent and from 14.3 to 30.8 percent 
for the "Circle Master" and the "Valley 1160", respectively. Where the 


the effect of wind was not included in the model the mean absolute de- 
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viations were increased to 17.5 percent for the former C.P. system 
as compared to 19.4 percent for the latter system. Close agreements 
were obtained between the measured weighted coefficients and those sim- 
ulated in the model (a maximum deviation of approximately 6 percent). 

Effect of wind on uniformity of water distribution was shown 
by the curves of the theoretical simulated application depths. The 
uniformity coefficients were affected by both the wind velocity and 
wind direction, however, the effect of wind direction on the uniformity 
was more pronounced. For a wind velocity of 20 mph and 0 degrees 
direction, where the line of measurement was directly opposite to the 
wind, the coefficients of uniformity were 87.3 percent and 91.4 percent 
for the "Circle Master" and the "Valley 1160" C.P. systems, respectively. 
For the same wind velocity of 20 mph and 80 degree direction, uniformity 
improved as the coefficients of uniformity were 96.02 percent and 96.54 
percent for the former and latter C.P. system, respectively. 

Performance of the C.P. systems was determined from field measured 
application depths as well as those simulated by the model. In general 
the field coefficients were lower than simulated coefficients, but both 


the field and the simulated coefficients were within the acceptable level 


of 80-90 percent. The spread of the weighted field coefficients varied 
from 87.6 percent to 93.6 percent for the "Circle Master" and from 88.4 
percent to 92.6 percent for the "Valley 1160" C.P. system. The weighted 
Simulated coefficients varied from 91.8 percent to 96.1 percent for the 
former C.P. system and from 93.4 percent to 96.1 percent for the latter 
C.P. system. The weighted coefficients, field and simulated, indicated 


higher uniformities than the respective non-weighted coefficients. 
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CHAPTER I 
INTRODUCTION 

Centre pivot (C.P.) irrigation is presently an increasingly popular 
method of irrigating large agricultural areas. The main advantages of 
C.P. systems are a high degree of automation, adaptability to various 
soil types and topography, and a higher application uniformity capa- 
bility. A relatively high initial investment is therefore justified 
due to the low labour requirements, the potential of increasing yields 
and high water application efficiencies. 

By definition, a C.P. system consists of a lateral line of sprink- 
lers continuously rotating about a pivot. Water under pressure is 
Supplied to the lateral at the pivot. The lateral is supported by a 
number of self-propelled towers, each having a driving mechanism utili- 
Zing wheels or tracks. The adjustable speed of travel of C.P. systems 
allows the application of 1/2 to 4 inches of water per revolution. Stan- 
dard C.P. systems are approximately 1/4 mile long, irrigating 135 acres 
of a 160-acre field. 

The important features of a well-designed centre pivot system, as 
recognized by the American Society of Agricultural Engineers (1), are: 
capacity to supply an adequate amount of water to the soil at the peak 
water use periods, to replace moisture in the soil profile frequently 
enough to maintain conditions for optimum plant growth, and to refill 
the soil moisture reservoir uniformly. 

The first two features are directly related to the soil type, the 
crop grown and the area to be irrigated. The usual approach to achieve 
relatively uniform water distribution is to increase the discharge from 


the first to the last sprinkler, that is, by increasing the application 
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rates in proportion to an increasing area irrigated by individual 
sprinklers. Previous research (4, 26) has indicated that the theoret- 
ical uniformity varies from the actual field uniformity. Large uniform 
ity deviations are the result of inadequate sprinkler systems design 
and external factors, mainly wind. The effect of wind on the distri- 
bution from stationary sprinklers systems has been documented (14). 
However, the distribution from C.P. systems as affected by wind has 


not been examined to any great extent. 
This study is an attempt to analyze experimentally and theoretically 


the interrelationships between wind, sprinkler patterns and uniformity of 
water distribution and to provide a model for predicting the application 
depths from centre pivot systems under windy conditions. 
Objectives 
The objectives of this study were: 
1. To evaluate the performance characteristics of two centre 
pivot sprinkler irrigation systems operating under normal 
climatic conditions. Field measurements of application 
depths and a continuous record of meterological data 
monitored were used to examine the effects of wind on the 
uniformity of water distribution. 
2. To obtain the desired wetted patterns of the sprinklers used, 
by simulating the flight path of water droplets emitted. 
3. To incorporate the predicted wetted patterns into a mathe- 
matical model developed to simulate the performance of the 
C.P. system operating under windy conditions. 
4. To validate the model by comparing the predicted and measured 


application rates. The effects of continuously varied wind 
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conditions on the uniformity of water distribution were 


examined using the model. 
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LITERATURE REVIEW 


The uniformity of water distribution by sprinklers has been of in- 
terest to agricultural engineers for the past half century. Research 
has shown that although sprinklers distributed water uniformly in a 
Circular area, the uniformity between the sprinklers was not satisfac- 
tory. The use of overlapping sprinkler pattems was an important break- 
through in solving the problem of poor uniformity encountered with cir- 
cular sprinkler water distribution. For stationary sprinkler systems, 
Christiansen (6) demonstrated that nearly uniform distribution is possible 
with proper sprinkler patterns and proper spacings of sprinklers. The 
development of highly automated sprinkler systems, such as centre pivots, 
further increased the potential application of sprinklers in irrigation. 
2.1 Measures of Uniformity 

Several expressions have been developed for the purpose of compar- 
ing the sprinkler distribution patterns. Christiansen ( 6) introduced 
the coefficient of uniformity (Cu) expressed as a percentage. The uni- 
formity coefficient (Cu) given below is dependent on the mean deviation 
of water distribution. 

2 (x; _ x) 
Cu = 100 ( LOS ares err Crelin’) 
Nx 
where 

Xj — individual readings 
R = the mean value of xj 
N - number of observations 


Wilcox and Swailes (26) derived a similar coefficient of uniform- 


ity but squared the mean deviations. The distribution, as expressed by 
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the following equation, was considered to be adequate for values 


greater than 70 percent. 
v= 100 (1.0 - a 


x 
where 


U - coefficient of uniformity 

SD - standard deviation of all observations 

x - the mean value of all observations 

Benami and Hore ( 2 ) proposed a new coefficient of uniformity 
by considering the deviations from the group means rather than the gen- 
eral mean. The general mean was computed from all observations. The 
group means were determined from observations above and below the gen- 
eral mean. This approach gave more weight to readings with extensive 


deviations from the general mean. 
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where 

Ao - uniformity coefficient 

Na - number of observations above the general mean 

N, - number of observations below the general mean 

T, - sum of the observations above the group mean in the above 
group 

T, - sum of the observations below the group mean in the below 
group 


D, - difference in the number of observations below and above the 
group mean in the above group 
Dp —- difference in the number of observations above and below the 


group mean in the below group 
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M, ~ group mean of the observations larger than general mean 

M. - group mean of the observations smaller than general mean 

Tezer (23), in investigating the reliability of the above three 
coefficients of uniformity, suggested a new coefficient (A), where 
(A) is the coefficient (A) expressed as a percentage (A = 0.6A,) - 
The coefficient (A) proved to be a better estimate of the degree of uni- 
formity than the coefficients (Cu, U) described earlier. This was due 
to the sensitivity of the coefficient (A) to variation in distribution, 
and thus gave a greater emphasis to extreme readings below the mean. 

A frequently used measure of sprinkler performance is the USDA 


pattern efficiency B.) (7), defined as the ratio of the average of the 


lowest 25 percent of the appiication depths to the average depth of 


application. 
(1/4)N 
Xi 
i 
aL 
Er, = 100 : 
(1/4)N x 
where 
sesece@. |OWest 25) percent of the readings 
xX - average of all readings 
N - total number of readings 


Skewness of the distribution either positive or negative which occurs in 
windy conditions would affect the above USDA pattern efficiency Epy, since 


it depends on the lower tail of the distribution. 


2.2 Factors Affecting Uniformity 
The main advantage of sprinkler irrigation is the capability to 
distribute water more uniformly onto the soil than any other method of 


irrigation. The factors affecting uniformity can be classified depending 
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on their origin as (19): 

a. Sprinkler head factors; represented by the size of sprinkler 
nozzle, type of nozzle, angle of nozzle, rate of sprinkler 
rotation and operating pressure. 

b. Distribution system factors; such as sprinkler and lateral 
Spacing, height of sprinkler above field surface, stability 
of sprinkler riser and variation of pressure along the lateral. 

c. Climatic factors; mainly wind speed and direction. 

d. Management factors; varies with type of sprinkler system. 

With stationary systems these are the system layout and dura- 
tion of system operation. With self-propelled systems the 
factors are the speed of the lateral and the alignment of the 
lateral and the sprinkler risers. 

Christiansen (5) made extensive field tests to investigate the 
factors affecting the water distribution of sprinkler systems. He inferred 
that the water distribution from sprinklers could be classified into 
Six characteristic geometrical patterns, designated as A, B, C, D, E, and 
F curves. ‘The A, B, and C curves varied from triangular to conical 
shape and gave uniform distributions when sprinklers were no further a- 
part than 55 to 60 percent of the diameter of a sprinkler wetted pattern. 
The D, E and F curves, however, resembled a rectangular cross-section 
and yield high uniformity coefficients with spacing up to 45 percent, 
dropping sharply between 45 to 75 percent and increasing again between 
75 and 80 percent of the diameter covered. 

Wilcox and MacDougald (25) selected eight types of distribution 
curves to determine their impact on the uniformity of water distribution 


and reported similar findings to those of Christiansen (5). The best 
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distribution curve was found to be one with a steady decrease in appli- 
cation rate from the sprinkler towards the boundary of the pattern. For 
Squared spacings the sprinklers should not be further apart than 60 
percent of the diameter of throw. A squared spacing results in better 
uniformity than a rectangular spacing for the same area covered when 


no wind is considered. 


Sprinkler wetted patterns are distorted when inadequate water 
pressure is supplied. Bilanski and Kidder ( 3 ), conducting laboratory 
experiments on medium pressure sprinklers, Showed that the uniformity 
of water distribution increased with rising pressures. With higher 
pressures greater breakage of the water stream occurs resulting in 
smaller droplets, filling the gaps normally caused by lower pressures. 
A further increase in uniformity was observed from sprinklers operating 
with oscillating arms. For a given operating pressure the distribution 
improved with a larger nozzle size. To minimize the decrease in uniformity 
due to variation in system pressure, Christiansen ( 6 ) recommended 
tnat the frictional losses along the stationary lateral should be kept 
Within 20 percent of the average operating pressure. 

The uniformity of water distribution is affected by wind be- 
cause of the skewed wetted patterns. At low wind velocity (less than 
3 mph) most sprinkler systems operate favourably. However, with a high 
Wind velocity (10 to 14 mph) sprinkler wetted patterns become asym- 
metrical, as a result of higher concentration of water near the sprink- 
ler and on the lee side. Christiansen ( 5 ) showed that with closer 
sprinkler spacings and adequate overlap between sprinklers the effect 
of wind on the uniformity of water distribution was less significant 


than that of other factors. Korven (14) in his testing various 
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revolving sprinklers, indicated that the effect of wind at high veloci- 
ties on uniformity of water distribution was quite dramatic. From the 
results, it was concluded that the average coefficient of uniformity of 
all sprinklers tested fell from &2 to approximately 32 percent when the 
wind velocity incresed from 4 to about 17 mph. At fairly high wind 

velocities the uniformity largely depends on the design characteristics 


of the individual sprinklers. 


2.3 Factors Affecting Sprinkler Drop Distribution 

The drop size distribution of sprinkler sprays and dynamics of 
the individual drops are of practical importance in sprinkler irriga- 
tion. Small droplets are subject to wind drift and evaporation loss- 
es which in turn result in distorted patterns, poor uniformity and 
reduced sprinkler efficiency. At impact, large drops with greater 
kinetic energy have the potential of breaking up the soil structure 


and also damaging the leaves of young plants. 


2.3.1 Sprinkler Jet Breakup 

A jet of water emerging from a nozzle under relatively high pres- 
sures breaks into small droplets because the state of equilibrium 
has been disturbed. 

Merrington and Richardson (15) summarized previous research con- 
cerned with different phases of disturbance of a liquid jet. In the 
first phase of disintegration, the jet becomes varicose due to the 
action of inertia and surface tension. The rate of increasing distur- 
bance is independent of the speed of efflux. In the second phase of 
disturbance, the jet is further deformed and the action of the air 


becomes more important than surface tension. Since the air resis- 
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tance increases rapidly with speed, the rate of jet disturbance 
increases with increasing speed of efflux. The third phase of jet 
breakup, in which the jet is completely disrupted, is termed "atom- 
ization" by Ohnesborge, as cited by Merrington and Richardson (15). 

Miesse (16) modified the criteria originated by Ohnesborge by 
determining the region where the three types of jet disintegration 
occurs. Based on dimensional analysis, Miesse indicated that the third 
phase of disintegration occurs in the secondary atomization region 
rather than in the third region. The prime forces acting in that 
region controlling breakup are viscous, inertia and capillary 
forces. 

Rouse et al. (20) described the jet breakup just after emergence 
from the nozzle as a result of continuous action of countless numbers 
of turbulent eddies no longer confined by a conduit. At this stage, 
the density of air is too low to have any significant effect upon 
a jet. The action of air begins only after the jet is sufficiently 
disrupted. The resistance of air is proportional to the square of the 
velocity or the jet and the cross-sectional area. The resistance 
rapidly decreases as the distance of travel of the eddies carrying 
water from the main stream forces them to spread laterally at dimin- 


ishing velocity. 


2.3.2 The Drop Size Distribution 

Merrington and Richardson (15) found experimentally that the mean 
drop diameter of water particles depends on the fluid viscosity and 
the jet velocity. The mean diameter was also found to be inversely 
proportional to the relative velocity between the jet and the surround- 


ing air. Considering the cross-sectional profile of the jet's relative 
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velocity to the air, the smallest relative velocities are near the 
jet's periphery. This will result in smaller particles at the peri- 
phery of the jet and the largest particles near the core. 

Seginer (21), examining the interrelation between the sprinkler 
drop size distribution, drop velocities and sprinkler distribution 
patter, made the following observations: 

a. With increased operating pressure, the mean drop diameter 

and the application rates decreased, the radius of water 
throw increased and the distribution pattern improved. 

b. Higher rates of sprinkler rotation tended to decrease the 

radius of throw. The rotation of the oscillating arm did 
not significantly affect the distribution characteristics. 

c. At a given distance, larger particles were formed from ro- 

tating sprinklers as compared with non-rotating sprinklers. 

Kohl (13 ) demonstrated experimentally the effect of pressure and 
nozzle sizes on the mean drop diameter. When the nozzle size was held 
constant and pressure varied, the increased pressure resulted in a greater 
number of smaller drops. Since jet velocity is proportional to the pres- 
sure, the experiment confirmed the relation of decreasing drop size with 
increasing relative velocity between the jet and the air. Similarly, by 
decreasing nozzle size at constant pressure the mean drop size is de- 
creased, but the effect of rising pressure on the mean diameter was more 
pronounced. Smaller sprinkler nozzles will produce smaller mean drop 
Sizes than large sprinklers under normal operating conditions. Operating 
small sprinkler nozzles at low pressures may result in larger mean drop 
diameters than the mean drop diameters produced by larger sprinklers at 


higher pressures. 
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Several methods of measuring droplet size have been employed to 
determine the distribution. The most common method of measuring the 
particle sizes is the "stain method" described by Hall (10). This 
method is based on the assumption that a drop falling upon a uniform 
absorbent surface leaves a stain of a certain diameter, proportional to 
the diameter of the drop. Seginer (21), using the stain method, de- 


termined the drop size from this relation: 


d = 0.444 p2/3 fon) 


d - diameter of the drop (mm) 


D - diameter of the stain (mm) 


2.3.3 Behaviour of Water Drops Emitted From a Sprinkler Nozzle 

After the initial disintegration of the main stream, the individual 
water droplets will deform or may further break up under the high velo- 
city depending on the size of the drop. Blanchard, cited by Umback and 
Lembke ( 24 ), concluded that drops below 4.6 fee quite stable, and 


drops between 4.6 and 5.4 mm break up upon shock. The trajectories of 


the drops will largely depend on their tangential velocities. The forces 


exerted on the particle are the resistance of the air, gravity and under 
windy conditions, the wind. 

Seginer ( 22 ) discussed some phenomena related to the tangential 
velocity of sprinkler drops. He assumed that the two forces acting on 
a drop are the drag force and gravity. The air drag force was described 
as an "ever increasing function of velocity (V) for any given body". 
The relation between the velocity and the drag force was expressed by 


means of the drag coefficient Cd. 
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where 

p- the air density (Kg - 52/m°) 

A — the projected area of the particle (m?) 

F'— the drag force (Kg) 

V - the relative velocity of the particle (m/s) 

The drag coefficient of a water drop was considered to be a func- 
tion of the Reynolds number and Weber number because of the deformation 


of drops at high velocities. At high velocities, the drag coefficient (Cd) 


approaches a minimum value and the drag force may be considered directly 
proportional to the square of the velocity. The effect of gravity on 


the particle depends on whether the drop has an upward or downward component 


of velocity. The drag force F was approximated as the lower function 


of the velocity V. 
= 4 
ce (4) 
where 


Cy - coefficient depending on value of n 


n - positive dimensionless constant 


The effect of the constant n on various trajectories of drops was 
investigated. A value of n = 2 proved satisfactory, on the average, in 


predicting sprinkler droplet paths. The equation of motion is then: 
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g - acceleration due to gravity (m/s°) 
+ —- sign indicating upward or downward movement 
Cc mua 
2 - drag coefficient 
2 - relative velocity of a particle (m/s) 
a2 
= acceleration of a particle (m/s?) 
dt 


Okamura (17) derived an equation of motion for water particles 
emitted from a sprinkler nozzle in still air conditions and under windy 
conditions. The drops were assumed to retain a spherical shape for the 
average flying time. Under no wind conditions, the equations of motion 
was based on the air drag force and the evaporation losses. The air 
drag force consists of two components; namely, the viscous and inertia 


resistance expressed respectively as: 


fea Oe eV a (6) 


1 2 


where 


7 - dynamic viscosity (Kg - s/m2) 

V - velocity of the particle (m/s) 

a - radius of sphere (m) 

Cd - air drag coefficient 

S - projected area of the sphere (m2) 


p - density of the fluid (kg - s2/m4) 


When the Reynolds number is greater than 1, the viscous resistance 
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is negligible and the inertia resistance may be considered the drag 
force. The drag coefficient Cd was approximated from the following 
equations: 


R,<100 C= 33.3/R, - 0.0033 R, + 12 
100<<R,<1000 ¢, = 72.2/R, - 0.0000556 R, + 0.48 p (7"4/brc) 


Ry = 1000 Ca = 0,45 


The equation of motion of the spherical drops in still air in the 


direction x and y are: 


Par3 
a Pa 3 . 


where 


X - acceleration in the x direction (m/s?) 


y - acceleration in the y direction (m/s2) 
Pa - density of the air (Kg-s?/m4) 

Pw - density of the fluid (Kg-s2/m4) 

D - diameter of the sphere (m) 


gq - velocity of the particle (m/s) 


The effect of evaporation causes a decrease in droplet size while 
the drop travels in the air. Okamura (17 ) showed that the change in 
diameter due to evaporation can be calculated from equations by incorp- 


orating the heat balance term in the equations of motion. Nevertheless, 


when particles greater than 0.5 mm are considered the effect of evapora- 
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tion may be neglected. 

The difficulties inherent in the study of the motion of particles 
when influenced by wind may be considered in two aspects; the lack of 
research on motion of water drops in still air and the problem of obtain- 
ing reliable data due to the turbulent flow nature of wind near the sur- 
face and the chaotic variation of wind in most climatic conditions. 

In deriving the equations of motion of sprinkler droplets under wind 
conditions, the following assumptions were made ( 18 ): 

a. Within some period of time, the wind velocity is assumed to 

be constant. The instantaneous values of the wind velocity 
are the sum of the mean and fluctuating values. 

b. The wind profile is assumed to be a uniform, vertical dis- 
tribution at a large distance from the surface. Near the 
ground the profile is described logarithmically, thus the 
Slope depends on the surface roughness. Coriolis forces 
have a negligible effect at this relatively small scale. 

Using the above assumptions, the equations of motion of flying drop- 

lets were derived for the case when the wind direction is parallel to 
the x axis of the x, y and z_ coordinates and the Magnitude of the mean 


wind speed is w. 
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Z = Ca p., 4D 


where 


ee , . . ; 2 
x - acceleration in x direction (n/s_) 
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z - velocity in z direction (m/s) 
w - wind velocity (m/s) 


g - acceleration due to gravity (m/s2) 


The trajectories of flying water droplets with various diameters 
emitted from a sprinkler nozzle are obtained by solving the differential 
equations (10), (11) and (12) for the known initial conditions, using 
the computer. From the calculated trajectories, drop size distribution 
at any level and at any height were evaluated. The phenomenon of 
wind drift was clarified by considering the average wind velocity and 
the varied discharge direction in respect to the wind direction in the 
equations of motion. Because of the presence of a vertical component 


of wind velocity, fine water droplets are carried away by wind drift, 


thus accounting for losses in sprinkler spray. 


2.4 Theoretical Distribution of Water by Moving Sprinklers 

In modelling the cross-sectional distribution of sprinklers, two 
types of distribution curves, resembling a triangle or half an ellipse, 
are commonly used. Both these geometrical curves were found to give good 
uniformity within specific spacings (4 ). The maximum rate with the 
triangular distribution is double the maximum application rate with the 


elliptical pattern. 


2.4.1 Single Moving Sprinkler 
Bittinger and Longenbaugh ( 4 ) analysed elliptical and triangular 
patterns of a single moving sprinkler to determine the degree of skewness 


when the sprinkler is moved in a straight line and in a circular path. From 
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the geometry, (see figures (]) and (la) , precipitation at any 


point (P) for triangular and elliptical pattern respectively, were 


defined as: 
eg 1 a (13) 
Dp =3 be = als (14) 
where 


h - maximum precipitation rate at middle of the pattern (in/hr) 
r - radius of sprinkler patter (ft) 


S - distance from sprinkler to point P (ft) 


The total depth of application from one sprinkler was evaluated by 
integrating the rate over the total time. The total depths from one 
sprinkler moving in a straight line at constant velocity assuming the 


triangular and elliptical patterns are respectively: 


2 the 2.1/2 _ 2 4 (i-m@)1/2 41 
Pex = Ba ( an Sp ars ae (15) 
ee eee Tea att) 
el 
2WR (16) 


Dy, ~ total application depth with a triangular distribution (in) 


De] - total application depth with an elliptical distribution (in) 


6p) 
I 


distance from the pivot to a point p_ (ft) 


Q 
I 


angle of rotation about the pivot (radians) 


@ - angular velocity of the sprinkler (radians/hr) 
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Figure 1. Application rate patterns of a stationary sprinkler. 


Figure la. Geometric Representation of a sprinkler moving in a 


circular path at constant velocity. 
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m-—- ratio ( 


R - distance from the pivot to the sprinkler (ft) 

@ - angle of integration; @ = k("%-a) radians 

T - time required for one half of the sprinkler pattern to pass 
DOMerE (hie) 

h - rate of precipitation at centre of sprinkler pattern (in/hr) 

r - radius of sprinkler patter (ft) 

The total depths from one sprinkler moving in a circular path 


assuming a triangular and elliptical pattern are respectively: 


GWT W/2 
2 Pueh 2 Deel 
D_ = 2ht os ff (s ne as 2 RS cos a | da en) 
0 
WT 
OF) 72 
4h - ney ff Ge i a 
Oe ere | ern, (1 rae sin 6) dd (18) 


Bittinger and Longenbaugh (4) observed the distribution for both 
pattems moving in a straight line by plotting (D, /2) and (D.4/%) against 
(m) from equations (15) and (16). For the case of a sprinkler moving in 
a circular path equation (17) was solved as an elliptical integral of 
the second kind E(k,%@). Equation (18) was not readily evaluated except 
by numerical means. They also noted that for (n) greater than 5 the cir- 
cular path could be approximated with a straight line path. 

The maximum precipation rate from a stationary sprinkler for both 
distribution patterns is near the centre. When the sprinkler is moved in 
a straight line or a circular path, this rate is further intensified. 

The pattern is skewed toward the centre of curvature when moved in an arc. 


The degree of skewness varies. inversely with the distance between the 
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sprinkler and the pivot. With effective overlap between adjacent sprink- 
ler and the pivot. With effective overlap between adjacent sprinklers 
both the triangular and elliptical distribution curve will result in 
high uniformities. The optimum spacings for the two types of distribu- 


tion curves respectively are 1.0 and 1.4 times the pattern radius (4). 


2.4.2 Rotating Lateral 


The application rate on a point (P) at a distance (S) from the 
centre pivot includes discharge from more than one sprinkler and can 
be determined by summing the individual rates from overlapping adjacent 
patterns. The equations for application rate for the triangular and el- 


liptical patterns, respectively are: 


eae: r 
fas 00) i, oath Say nah a 
Ae =)> hy = (S° +R. - 2R,S cos @) (19) 
i i 
we = > HE ae ~ 5? - Re + 2R.S cos a)? VA 
el nee at i ut 
o a 
where 
i - subscript referring to gen sprinkler on the lateral 
A, - application rate with a triangular distribution (in./hr) 
A_ - application rate with an elliptical distribution (in./hr) 


e 


h,r,S, R, , - aS previously defined 


The equations (19) and (20) solved by Heermann and Hein (11) for small 

increments of (@) from zero to (WT) gave one half of the symmetric 

rate-time relationship for a complete pass of the patterm over point (P). 
For the C.P. system analyzed, the equations solved to determine the 


total depth at a distance (S) from the pivot for the triangular and el- 
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lliptical patterns, respectively, were: 


= h T 
- fs 4 
Oh 2. h,T, 2 = | e (e + Re =) FARES Ve O'S: a da@ (21) 
i=1 se a 
0 
d 
N WT. 
2 4n. (n.+ m.) 25 
1. 
Pa =e Lana fo (1-2 ain’o) ae 
1i=1 0 tec 


(22) 


The distribution patterns for a C.P. system are obtained by nu- 
merically solving equations (21) and (22). The large number of re- 
petitive solutions suggests the use of a computer. For application 
depths beyond 650 ft (i.e., the skewness for the circular path is less 
than 0.01 in. from inner to outer radius for one sprinkler) straight 
line approximations for triangular and elliptical patterns could be 
used (11). 

For comparison of theoretical and field distribution Heermann and 
Hein (11) found a close fit between computed and measured application 
depths for triangular and elliptical distribution were almost identical. 
The small differences of the two distributions were associated with the 


close spacings of sprinklers in the systems studied. 
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CHAPTER Iii 


MODEL DEVELOPMENT 


The model simulating the centre pivot irrigation system consists 
of two separate parts. The first part involves the dynamics of water 
droplets emerging from the rotating sprinklers. The water droplet paths 
define the basic wetted patterns of a sprinkler for a range of angles 
and discrete values of particle sizes, initial particle velocities and 


wind velocities. The second part of the model, which simulates the rel- 


ative performance of the C.P. system, performs two distinct functions, 
namely: 

a. The interpolation of the geometric variables that define 
the finite number of calculated wetted patterns. Predicted 
wetted patterns could then be determined for a range of 
initial conditions encountered. 

b. The determination of the system performance described by the 


application rates, depths and coefficients of uniformity. 


3.1 Dynamics of Water Droplet in Flight 


Experimental research on the distribution of water droplets inside 
the wetted patterns has shown that the largest droplets assume the 
largest flight trajectories (21). Blanchard, cited by Umback and Lembke 
(24), indicated that, after the initial break-up of the main stream, the 
water droplets will further deform or disintegrate depending on the tan- 
gential flight velocities and the droplet sizes. In modelling the 

the flight of water droplets from a sprinkler the following assump- 


tions were made: 
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a. The viscous component of the drag force acting on the water 
droplet is negligible since the Reynolds number under most 
flight conditions is greater than 1. The drag coefficient Cd 
is therefore a function of Reynolds number. 

b. .The water droplet retains a spherical shape throughout 
the flight time. 

c. The effect of evaporation on droplet size is negligible. 

d. For the flight time considered, the wind velocity is con- 
stant with a uniform vertical profile. The wind resistance 
acts only in the horizontal direction. 

Taking the above assumptions into consideration, the equations 
of motion of flying water droplets were derived. The forces acting 
on a water particle at a particular instant in flight are the drag 
force (Fa)s the gravitational force (mg), and the inertia force (mx) 
as shown in figure 2, Applying D'Alembert's principle, the sum of 
forces on the particle in the inertial reference frame (R) equals 


Zero. 


Figure 2. Forces acting on a water droplet in free flight. 
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The drag force acts in a direction opposite to that of the relative 


velocity between the particle and the air and is given by the 


equation: 

p= “KW 
and 

k = 1/2 Ca pA 
where 


C a the drag coefficient (dimensionless) 


P - the density of air (slugs/ft>) 


A - projected area of the spherical droplet ( £t*) 


Consider a water particle at an instant (t) in flight whose 
position relative to the ground was described by a coordinate system 
R(x,y,z) aS Shown in figure 1. Assume x, y and z to be the positive 
particle velocity components in the positive x, y and z directions. 
The velocities of the water particle (p) and wind (w) were given as: 


VS xe yin eek 
Oba “Rr Pp 
W 


¢ e 
Xx it j 
WwW Yy,J 


The relative velocity between the water particle and the wind is: 


Eee a ee a 
el SN ee 2 aku Day al / 2 
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witha, and Y defined as the angles the relative velocity vector 
Vp makes with the x, y and z coordinates respectively; the drag 


force can be resolved in x, y and z directions to obtain: 


a aoe =< 4 i 
F d fi + f fix 
f = Ky" cos @ 
pe R 
t = -kVS COs B 
ie 2 
f. = kV, cos ¥ 
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where the directional cosines are: 


cos @= s - of VR) 
cos B = ea eee Ae, 
cos y= (2/V py 


Substituting fora, and y the components of the drag force are: 


a = “KV, - x ) (7.4) 
f, = “KVR (Yp ae) (8%) 
cS = “kV, &,) (795) 


Summing the forces acting on the particle we obtain the following 


reactions: 


F, = 0 = -mz_ - kV,z -mg 


The accelerations of the particle in the x, y, z direction are: 
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The differential equations (10, 11, 12) were solved for known 
initial conditions by means of a continuous simulation modelling 
program (CSMP). The model performs continuous integration of the 
particle acceleration and velocity during the specified flight 
time (t ). The simulation model not only describes the trajectory 
of the particle in flight but also permits the geometry of the wet- 
ted pattern to be defined when the water particle strikes the ground 
surface. 

The ratio (k/m) must be known to solve the differential equations. 


The mass of the particle (m) for a predetermined particle diameter 


(d) is obtained by using the earlier assumption that the droplet is 


spherical in shape. Therefore, 


2 
A T (13) 
and 
ae AG) 
Ts g (14) 
where, 
A - projected area of the spherical droplet (£7) 
d - diameter of droplet (ft) 
g - acceleration due to gravity GY) 
ve - volume occupied by spherical droplet (£3) = = 1a 
fe - density of water (62.4 1b/£t”) 
G, - specific gravity of water (dimensionless) 
The mass density of the air (fP) is obtained from the relation: 
p = = (slugs/£t?) (15) 
where, 
Vv - specific volume of air at 65°F 
The drag coefficients (C5) were estimated from the diagram 
as a function of Reynolds number Re: 
R= & (16) 


where, 
v - kinematic viscosity of air for atmospheric pressure 
fry ae radretarper attire Or GS n(re aia) 
Various values of the ratio (- < are determined for selected 
range of discrete values of particle diameters (d) and the associ- 


ated initial velocities by using the above relationships. 


Initial analysis indicated that the simulated patterms could 
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be categorized into two types: 

a. when no wind condition existed, the wetted patterns ap- 
proximated a circular shape that was described by the 
equation of a circle Ge + yo = r’). 

b. with the inclusion of wind as a vector in the initial 
conditions, the resulting skewed wetted patterns were 
best approximated by joining two halves of two ellipses 
with a common minor axis og fae + y*/b- = 1 and 
jer + yo Agim 

These geometric variables (a, b, c, e, and x) were simulated 


in the model to represent all possible wetted patterns that resulted 


from the initial conditions considered in this study. 


3.2 Determination of Predicted Wetted Patterns by Interpolation. 
The .basic wetted patterns that were simulated in the first part 
of the model were defined by the variables a, b, c, e and x; where a, 
b and e are the semi-axes of the elliptical shaped patterm, c the loca- 
tion of maximum application rate relative to the sprinkler and x is 
the radius of the circular pattern when no wind condition prevail. 
(See fig. 3). For specific values of the system variables such as 
sprinkler nozzle angles (Y), sprinkler height (z) and wind velocity 
he these geometric variables (a,b, c, e and x) were interpolated 
from those simulated with the measured particle size (P.) and 
initial particle velocity (Vea! as the specified arguments. A two 
dimensional interpolation technique ( the natural bicubic spline 
approximation) was used, in this section of the model, to interpol- 
ate the given set of unequal spaced data (9). A subroutine, ICS2VU, 


from the IMSL statistical package (12), was used to do the interpol- 
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ation. The subroutine was necessarily incorporated in the model to 
eliminate the lengthy computer time required in simulating the predicted 
wetted patterns for all combinations of initial conditions. 

Preliminary investigation of simulated variables a,b, c, e and 
xX Suggest a linear relationship for a change in wind velocity. There- 
fore, the predicted wetted patterns for measured wind velocity 
within the simulated range were approximated by linear interpolation. 
The predicted wetted patterns outside the simulated range of wind 
velocities were extrapolated assuming that the linear function 


holds. 


NO WIND 


SPRINKLER 
LOCATION (Q) 


<a WIND 


Figure 3. Geometric variables and configurations describing the 
sprinkler patterns. 
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3.3 Performance Parameters of C.P. Systems 

This section of the model evaluates performance parameters 
of a C.P. system. The relative performance is assessed in terms of 
application depth (D.), application rates (A_) and the uniformity 
coefficient (C)- In order to determine these performance parameters, 
the dimensions (a,b,c,e and x) of the predicted wetted patterns 
and the system characteristics were used in the program proper. The sys- 
tem characteristics were the maximum application rates for the triangu- 
lar distribution (H,) and gee elliptical distribution (H3)i, pene Sprinkler 


c 


volume discharge (V) and the angular velocity of the centre pivot (@). 
In modelling the distribution of individual sprinklers along 
the lateral, triangular and semi-elliptical distribution profiles 
were assumed from earlier research (4). The two distribution 
profiles differ in that the maximum application rate of the ellip- 
tical case is approximately half that of the triangular profile. 
The location of maximum application rates for circular wetted patterns 
is located immediately below the sprinkler, but in the case of skewed 
wetted patterns, however, this is offset a distance (C) from the 
sprinkler (Q) along the major axis. The assumed distribution profiles 
for both circular (no wind condition) and skewed wetted patterns 


are shown in figure 4. The maximum application rates for circular 


wetted patterns (no wind condition) are given by the following: 


y= 289V 
eta ye a7) 
aoe 45.94V 
e 2 (eS) 
x 
where: 
H, - maximum application rate of circular wetted pattern with 
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triangular distribution (in./hr) 


ae - maximum application rate of circular wetted pattem with 


elliptical distribution (in./hr) 
V - discharge volume of sprinkler (gpm) 


x - radius of circular wetted patterns (ft) 
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Figure 4. Triangular and Elliptical distribution profiles under 
wind and no wind conditions. 
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For wetted patterns distorted by the effect of the wind, the maxi- 


mum application rates are: 
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Aw - maximum application rate of distorted wetted patterns 


with triangular distribution (in./hr) 


ee - maximum application of distorted wetted pattern with 
elliptical distribution (in./hr) 

a,b,e - dimensions of skewed wetted patterns (ft) 

In determining the application rates (Ap) for any given point 
Within the wetted pattern, the location of point P and the point of 
maximum application rate thy need to be described with reference to 
the centre pivot (0). The analysis considers the two following cases, 
depending on whether or not the distorting effect of the wind is included. 
Case I: Circular wetted patterns (no wind condition). 

Consider an instant (t) when the lateral (OM) makes an angle 
(@) with the (x) coordinate, and a sprinkler located at point (Q) 


along the lateral imposes a circular wetted pattern of radius(x), as 


shown in figure 5 . Let P be any point within the wetted pattern 


(T), and the distance (OP) in relation to the point (0) is: 


Qe =(R? + 5? - 2Rs cos (9-8) ) * C21 a) 


where: 


QP - distance of point P to the sprinkler location Q (ft) 


R - distance of sprinkler Q to the pivot O (ft.) 
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Dee =rcsuance: Of point P to the pivot..0 > (fe) 

§@ - angle between the lateral OM and the x coordinates (radians) 
P - point represents a measuring recepticle along a row of cans 
B - angle between OP and the x coordinate (radians) 


Figure 5. Circular wetted pattern of a moving sprinkler under 
no wind condition. 


The application rate of the sprinkler (Q) at the point (P) for both 


triangular and elliptical distribution are: 
H 
pa 


and 
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In simulating the rotating lateral the effect of overlap by adjacent 


sprinkler patterns is incorporated as: 


: Ae (i) 
eS Js OP 
tr ) Hy (i) TR mits (24) 
(1) 
ea 
n L 
H : 2 
2 are: 
A >; = Oe fh = wore ( 25) 
el (i) ae) (i) 
i=l 
where: 
1 = subscript referring to the pe sprinkler on the lateral. 
n = number of sprinklers with overlapping wetted patterns. 


The application rate at the point (P) varies continuously as the moving 

wetted patter passes over it in a given period of time. For a small time 
interval (A t) the application rate is assumed constant, the application 
depth at point (P) is hence a sum of the partial instantaneous application 


depths: 
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and 
n Ak 
172 
pote. ) a ( eee OESS pat (27 ) 
(i) (1) (i) (j 
i=l j=l 
where 
J - subscript referring to “ee time interval t. 


T - total time when point P is located within the wetted pattem (hrs) 


Case II: Skewed Wetted Patterns 


In this case the skewed wetted patter assumes an elliptical shape 
with wind acting in a direction parallel to the principal axis of the ellipse 
making an angle ( @ ) with the lateral (OM) (see fig. 6). The point of maxi- 
mum application rate (A) is displaced a distance (c) from the sprinkler lo- 
cation (Q) on the lateral. To compute the application rate of any point 
(P) located within the skewed wetted pattern, the dimension AP and AZ in 


relation to the pivot O are required. 
2S (v2 + v2 } ae (28 ) 


where: 


AP distance between the point P and the point of maximum applica- 
PionataArerA. sm(Et..) 
Uo = Sicos (6—Y) —R cos (¥+9 ).: (£b) 
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R  - distance between the sprinkler at Q and the pivot O. (ft) 


S  - distance between the point P and the pivot O. (ft) 
c - distance between point A and the sprinkler Q. (ft) 
6 - angle the lateral makes with the x coordinate of the pivot 
OF CEC) 
= (@+@). (radians) 


angle between the lateral and the direction of wind. (radians) 


angle between OP and the x coordinate of the pivot O. (radians) 


and, 
2a 
Seinke ay 1s) 
ae (3 5 — yee ( 29) 
a’ cos°-¢@ +b sin’¢ 
where: 
=< U 


a,b - the semi-axes of the ellipsoidal wetted pattern. 


The major axis (a) is equal to e in Equation (29) when the point of 


interest (P) lies above the minor axis (b). 


The application rates at the point (P) of the sprinkler at (Q) and 


from the adjacent contributing sprinklers are: 
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where: 


Arp? Agr, - the application rates at a point (P) from the pivot (0) 
for the triangular and elliptical distribution profiles, 
respectively. (in./hr). 

i - subscript referring to it sprinkler on the nee 
n - number of sprinklers with overlapping skewed wetted 


patterns at point (P). 


The total application depths at the point (P) for the cummulative time 


interval ( At) are: 


: sano, eae 
Do iz Du q pe) i" AZ,: ES A ce (228) 
i=] ee j=l 
n Tey 
, Ht (i) ( =? Ps, ) 172 
De = ye = Ai Nae Ne (330) 
peace: i= 
where: 


Drpr Der, ~ the respective application depths at point (P) for tri- 
angular and elliptical distribution profiles within a 
time span (T). (in.) 

T - total time when the point (P) is located within the 


wetted patterns 


The uniformity of distribution for the C.P. system was evaluated by 


means of the following coefficients of uniformity: 
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oe D,..- 
Cu ="100 int of AL aE te BG) ) (34 ) 


Cuyy) = 100 fl - ey i a | ) (35) 


where: 


Cu(q)s Cu (2) - the respective weighted and non-weighted Christiansen co- 
efficients of uniformity (percent). 
Di; - the application depth at point (i). (in.) 
- the distance from the pivot O to the point (i). (ft) 
i - subscript referring to the ith point along the line of 
measuring cans. 


N - the number of observations. 


These two coefficients of uniformity were incorporated in the model to 
assess their sensitivity in relation to the effects of variable windy condi- 
tions affecting the C.P. system. 

The application depths at each predetermined point along a given radial 
line were computed for both triangular and elliptical distributions. The 
system analysed, in this study, took into account the overlap effect of ad- 
jacent sprinklers along the lateral. The total depths at a distance (S) 


from the pivot were determined by solving the equations (30), (31) and (2), 
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(33) for one pass of the lateral at half degree increments (Aad= 3°). 
The overall application depth distributions for both distribution 
profiles are determined by using a numerical method in solving the 
equations. Considerable computer time is needed since a large number 
of repetitive solutions are required. 

The Fortran IV Computer Language was used to program the model- 
ling formulation simulating the dynamic operation of two manufactured 
centre pivot irrigation systems (Circle Master, Valley 1160) used 


in Southern Alberta. 
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Chapter IV 


CENTRE PIVOT SYSTEM SIMULATION MODEL PROCEDURE 


Sprinklers of the C.P. systems were modeled in the program using a 
CSMP technique to determine their respective wetted patterns. The model 
took into consideration the effect of wind. Predicted wetted patterns 
were then modeled in the main program to assess the performance of the 


functional C.P. system. 


4.1 Model Input Requirements 

The sprinkler simulation model was executed on the IBM 1800 digital 
computer using the CSMP package. ‘The program proper, however, was written 
in Fortran IV for use with the University of Alberta, Michigan Terminal 
System (MTS). Program execution was performed on the Amdal 470 computer. 
Once the model was programmed, the rapid speed of computation permitted 
the study of variations attributed to particularly sensitive input para- 
Tevers-: 

The information required in the sprinkler simulation model deals with 
the geometric and physical characteristics of the sprinklers and the in- 
fluencing environmental parameters. These parameters are listed as follows 
along with their corresponding variable names: 

i. Angle of sprinkler range nozzle (ALFA) 

2. Height of sprinkler above ground surface (ZH) 

3. Projected area of spherical water particle (AP) 

4, Mass of water particle (MAS) 

5. Initial particle velocity (IV) 


6. Drag coefficient (DCF) 


4l 


_s 
&: art by Sue my ; 


. fs Shed 1913 ca | fearete 
ave Sasa ee Whe ay St! pl us 


- i 
Natta S ¢ - - 
y 5 rs * 4h FF " red 


7 vo Ue B \ RIT (iW Ligdige en ees 


ie Leiicae: ‘ot te ‘n 5 al > EeyEti ae 


= 


j iy ‘ ee Oe | separ 
ed 7 


r 
s a ur Ste 
, ee : 
| E 3 
Gy 
{ C 
i Go Ue ~ris > e : aa 
~ ah) - 
= 
“ah - a fy a 
* eh bho * : 


{eas RSE) lls Sentara tx os 


| dah wate cewhei : 
. . (vr ae 


7. Density of air (DA) 

8. Acceleration due to gravity (AG) 

9. Wind velocity and direction (WV) and (WD) 

10. Total integration time and the desired integration interval (DT) 
and (TI) respectively. 

For the range of initial conditions considered in this study, the 
corresponding geometric variables (a, b, c, e and x) obtained as output 
from the above simulation program are stored as arrays in disc files for 
specific input ranges of particle size (xx) and particle initial velocity 
(yy). These geometric variables (a, b, c, e and x) were assigned variables 
TA, TB, TC, TE and TX in the storage files. 

The input information required in the main simulation program deals with 
not only the measurable environmental and particle trajectory parameters but 
also the system characteristics. The measured field parameters are: 

1. Particle diameters (PS) 

2. Wind velocities (W) 

3m Initial particle velocities (VIN). 

The program at this stage performs the necessary two stage interpolation 
to determine the predicted wetted patterns of the individual sprinklers. 

The system characteristics are next entered as: 

1. Location of sprinklers on the lateral (R) 

2. Number of sprinklers on the lateral (M) 

3. Sprinkler volume discharged (VV) 

4. Angular velocity of lateral (T) 

5. Direction of wind velocity relative to the lateral (AL) 

6. Location of measuring receptacles along a line relative to 


Pivot (Y) 
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7. Number of measuring receptacles along the line (MO) 

8. The angle between the lateral at any position and the x axis of 
a co-ordinate system established at the pivot O (THE) 

9. Integer specifying the number of increments of (A @ ) for constant 
wind velocity and direction (INC) 


10. Angle between the line of receptacle (OP) and the x axis (BETA) 


The computer model then completes the computation to determine the out- 
puts for C.P. system performance analysis. Several options are available 
in the model to determine the sensitivity of the input variables in an 


attempt to obtain optimum performance of the C.P. system. 


4.2 Model Operational Logic 
The sets of simultaneous differential equations, that describe the 
flight dynamics of a water particle emitted from a sprinkler, were solved 
by means of the CSMP model using the block diagram approach. The schematic 
diagram showing the functional blocks for this model is shown in Appendix 
(I). The program uses basically two Fortran statements, namely the configu- 
ration statements which define the inter-connection and function operation 
of the blocks and the parameter statements that associate numerical values 
with the block elements to particularize their functions. 
The C.P. system simulation model performs the following functions at 
various stages of the program: 
(a) Firstly, the program compiles the input data (i.e. the output 
information of the CSMP program, measured field parameters and 
system characteristics) from the source files using several read 


statements. 
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Interpolation is then used to determine the predicted 

wetted patterns for each individual sprinkler along the lateral. 
The predicted wetted patterns are defined for a range of measured 
field parameters and characteristics of the C.P. system consid- 
ered. When wind effects were taken into account, both two- 
dimensional and linear interpolations were performed; however, 

for no wind situations only two-dimensional interpolation was 
required. 

The model then proceeds to perform a search to determine whether 
the point (P}) is enclosed within the boundaries of wetted patterns, 
and the partial application rates for poth triangular and ellip- 
tical distribution (TR, EL) are calculated; then separately summed 
to compute the total application rates. The application depths 
for both distributions are computed as the product of the total 
application rates and the differential time ( At) required for 
the lateral to sweep across an angle of half a degree. However, 
if the point (P}1) is not enclosed by any wetted patterns the 
program continues to select the next adjacent (Pj) along the line 
(OML) , as shown in figures 5 and 6. A similar test and conm- 
putation is imposed on point P5, nevertheless, when five consecu- 


tive points (Pi, Pos ...- P.) are found to be outside the bound- 


5) 
aries of the wetted patterns, the search and test is extended 


to the points along the next lines (OM. ) and (OM, 


analysis was performed for 180 half a degree increments of 


ths 


the lateral position to ensure that the initial and 


final influence of the moving wetted patterms on the points 


along line (OM,), (OM) and (OM3) were fully considered. The 


i _ ie) 


a - i uaead ise 


wim 


ee 77) 1 & 
ae ce A ie 
% A elite Win eS aT 
E LY aa 1 ‘ & 7 
bis Nien Se? oe RN h oer yy:  ) res tcp > JE BEE 
I { 
sve see ehhh 
oO eee awe. cotter tion fea eo Gee os 
: ; _ wl a af 
7 
uhhh Sed ee. Sed Age 448 eal Aes: iy 
; : = Sih 
; ; a 
‘ ee 2 COM 
a : 
ne . a Oo rental At 
r a, | 4 “) af AEt4 abi. 
J 
oo 4 
i, Areracts hy wei rd 7 eect PLY - 
rie ‘> ie ; 
io 7 > Fs 
+3 itd 4 va Wie i ee _ Do 
an 
c - 74 « ‘ y! 
- A 
? ‘ bts SEF \ 4 ' ‘y 
ae hl \ ’ | 
, efi is ? 4 4 a 
’ 
; ul a), 7 { ef oie BS) Os oie: 
a 
4 "i tit 7 4 ry tad Seaues\ ee 
; ; . - er 
ear) |i , a nije ihe el at ij y ‘a = 
.7 - 
Toit ae Cit) arse tAg Hossam io byt a sil tesened gh if 
yy Aor 
iv r r ‘fn J 
74 fire 


° a) 
re nary “Ss 


total application depth for each point (P)) is computed as the 
sum of each depth calculated at every increment of the lateral 
position considered. 

(a) Finally, the coefficients of uniformity for the C.P. system are 


computed. 


The flow chart of the program simulating the C.P. system is shown in 
Appendix (II). After completing the computation of performance parameters 


for all desired input variations the program may be started over at the 


beginning to model a different set of system characteristics or terminated. 
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CHAPTER V 


COLLECTION AND ANALYSIS OF FIELD DATA 


Two different centre pivot sprinkler irrigation systems were selected 
in order that their relative performance could be evaluated. Both C.P. 
systems were located on sites that were free of wind barriers so that the 
micro-climate represented the typical climatic conditions of Souther 
Alberta. Perennial forage crops with rapid growth rates such as alfalfa 
and grass were grown on these sites. The frequent harvest of these crops 
permitted repeat measurement of application depths without the interference 


of the crops. 


5.1 C.P. System Characteristics 

The two C.P. sprinkler systems under investigation were 'Valley 1160' 
manufactured by Valmont Industries Inc.,Valley, Nebraska, U.S.A., and 'Circle 
Master' manufactured by Pierce Corp., Eugene, Oregon, U.S.A.. 

The two systems were selected because of their conceptual differences 
in design such as height of lateral, rotational speed of travel, sprinkler 
spacing, number of sprinklers and the type of sprinklers used. The general 
specifications of the systems are given in table l. The detailed per- 
formance data for the sprinklers used in both systems is given in Appendix 
( 3). The majority of these system parameters for either systems were 
fixed as specified by their manufacturers. Evaluation of variation in 
system performance attributed to change in system parameters were not 
possible since field measurements were taken on sites at the disposal of 
the owners. The rotational speed of system travel, however, is the only 


parameter adjustable to some degree. Two rotation speeds for each system 


were used to increase or decrease the application rates of the sprinklers. 
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5.2 Experimental Layout for Field Test 

The application depths for the C.P. systems are traditionally measured 
using a grid system of receiving receptacles. Conventional cans ( 4 in. 
in diameter and 6 in. in height), held in place by a 9 in. long rod were 
used to capture the water droplets. In each row the cans were spaced 20 
feet apart, the first can, however, was spaced 10 feet from the centre 
pivot. The accuracy of the measured application depths was enhanced 
by placing two other rows of cans set at an angle of 3 degrees apart. In 
an earlier analysis (section 3-3) the rows of measuring cans represented 
lines (OM,), (OMp), and (OM3), while the individual cans were designated 


as points Pj, Po... P Additional rows of measuring cans were placed in 


ne 
a Similar fashion at intervals of 90 or 120 degrees between the grids of 
cans as shown in figures 7a and 7b respectively. 

From meterological records the prevailing wind directions in 
Southern Alberta are west and south west. The grids of measuring cans on 
both sites were placed in directions as shown in figures 7a and 7b 
in order that the varied effect of the wind direction relative to the lat- 
eral could be readily investigated. 

Meteorological data was recorded at the station installed on the site 
sufficiently distanced away from the effect of the sprinklers. A thermo- 
hygrograph and a wind recorder were located in a standard meterological 
booth 4 feet above the ground. Temperature and relative humidity were 
measured by the thermo-hygrograph and recorded on a chart calibrated fora 
period of 7 days. Wind velocity and direction were sensed 9 feet above the 
ground by a weather vane and separately recorded on charts of an event 


recorder. The instruments functioned continuously during the experimental 


period. 
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Figure 7a. Schematic of experimental field layout at site l. 
(Cranford, Alberta. ) 
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Figure 7b. Schematic of experimental field layout at site 2. 
(Granum, Alberta.) 
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The amount of precipitation and evaporation during the experimental 
period were measured using a set of control receptacles. These receptacles 
were identical to the measuring cans used in application depths measure- 
ments and were located away from the effect of the sprinklers. In addition 
to this equipment, two tipping bucket rain gauge recorders were used in the 
experiment. The rain gauges were calibrated to respond to a change of 1/100 
of an inch of precipitation depth. Continuous recording for 24 hours was 
recorded on the charts. The tipping bucket rain gauge recorders were used 
to spot check the rainfall precipitation and to synchronize the measurement 


of application depths to the meteorological data recorded. 


5.3 Application Depth Measurements 

The measurements of application rates for both C.P. systems were made 
for two rotational speeds of travel during periods of varied weather condi- 
tions. After the lateral made a complete pass over the grid of cans, the 
volume of water collected in each can was measured using a graduated cylinder 
and converted to a depth in inches. When application depths measurements 
were delayed, the necessary corrections were made for either evapora- 
tion losses or gains due to precipitation. The loss or gain was estimated 
from the change in volume in the control measuring cans. For the hydraulic 
driven C.P. system (Valley 1160), some of the measuring cans were affected 
by the water discharge from the piston. Measurements in this case were 
discarded and an average reading was taken by considering the measurements 


of the adjacent cans. 


5.4 Sprinkler Droplet Size Determination 
The size of the droplets with the longest trajectories was determined 


for the purpose of simulating the sprinkler's wetted patterns. Droplet 
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sizes were measured for each individual sprinkler of both C.P. systems with 
the exception of the large volume 'end gun' sprinklers. In the case of 
sprinklers with double nozzles, droplets emitted by the range nozzle only 
were considered. The size of droplets were determined experimentally using 
the stain method (10). 

Filter papers (Whatman No. 1) sprayed with a suspension of methylene 
blue in carbon tetrachloride gave permanent stains when a water droplet fell 
on the surface. The diameter of stain was assumed proportional to the as- 
sumed spherical flight diameter of the particle. Calibration of equation (2) 
(section 2.3.1) for the filter paper, conducted in the laboratory with known 
diameters of water droplets, indicated (see figure 8) that the equation can 
be used in approximating the droplet size measured in the field. 

For field measurements of the emitted droplets a box containing the 
stain-sensitive filter paper was used. The box had a sliding surface that 
operated like a camera shutter. The sprinklers were held stationary when 
the filter paper was exposed to the farthest falling droplets. Repeated 
sampling was necessary for each sprinkler to ensure that the falling droplets 
did not give large overlapped stains. From these resulting stains the aver- 


age diameter of the droplets for each sprinkler was estimated. 


5.5 Field Data Analysis 

The variations in wind velocities and wind directions were recorded 
continuously throughout the experiment. Since it was not possible to deter- 
mine the application rates measured at each individual can in the row during 
a complete pass of the lateral, the effects of wind on each sprinkler cannot 
be fully evaluated. Therefore, a mean wind velocity and direction was used 


to relate the wind effects to the uniformity distribution of the sprinklers. 
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This mean velocity and direction was derived as an average of the values 
recorded for every 15 minute time interval throughout the duration of a 
Single lateral pass. 

Wind velocities and directions were required as input variables in the 
C.P. performance model as discussed in Chapter IV, (section 4.1). Because of 
the lengthy time taken for the lateral to make a complete pass across the 
grid of measuring cans (approximately 14 hours, 1/4 revolution), simula- 
tion of each change of the continuously fluctuating wind conditions was 
not feasible. For this reason, an approximate time interval of 15 
minutes was specified during which a representative average wind velocity 
and direction was computed for each interval throughout the experimental 
period. The time interval was adjusted to ensure an integer number of half- 
degree increments of the lateral rotation. This was necessary to facilitate 
the iterative computation in the main program. 

A range of particle sizes and initial velocities of the droplets emitted 
by the sprinklers was required as input variables in the simulation of pre- 
dicted wetted patterns. Sufficiently large samples of stains were obtained 
from field tests for all sprinklers to permit the determination of a range 
of particle sizes. The average sizes measured in the field were, therefore, 
assumed to always fall within this range. The range of initial velocities 
was determined from the manufacturer's specifications of the individual sprin- 
klers as shown in Appendix (IIIA) and (IIIB). The initial velocities were 


noted for a given discharge and nozzle size. 
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CHAPTER VI 


RESULTS AND DISCUSSION 


Simulation of the relative performance of the C.P. system was made 
possible by first determining the wetted patterns that were developed in 
the first part of the model ( see section 3.1, Chapter 3). The flight of 
water droplets was simulated for known initial conditions for each 


sprinkler along the lateral. 


6.1 Input Variables for Model Simulating the Flight of Water Droplets 

The initial particle velocities (IV) were calculated using the man- 
ufacturer's specifications (see Appendix III) as a guide. The droplet 
flight characteristics were determined from field measurements and applied 
in the theory that described the forces developed in the moving fluids. 
Preliminary estimation indicated that approximately 500 computing hours 
were required on the IBM II80 Computor to simulate the requisite trajector- 
ies. To reduce the large number of computations yet maintaining 
sufficient accuracy, the simulations were performed for all combinations 
of a particular range of input variables. This range of input variables 
enveloped all the situations that were encountered in the field. The 
interval of the discrete values of those variables were selected to 
maintain the accurracy. The droplet sizes, for example, were input as 
2,3,4, and 5 mm in diameter as shown in table (2), for initial droplet 
velocities of 90 ft/sec and wind velocity of 10 mph. The range of drop- 
let initial velocities (75,80,85,90,95,100,105 ft/sec) and wind velocity 
(0,10,20 mph)were used in the simulation model for given physical char- 
acteristics of each sprinkler. The resulting parameters defined the 


basic wetted patterns that were used to derive the actual wetted patterns 
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by method of interpolation. 
6.2 Wetted Patterns 


From the output of the simulation model the boundaries of the wetted 
patterns are geometrically described either by the equations of a circle 
of radius (x) (no wind condition) or that of two ellipses with a common 
minor axis. The validity of using these equations used to curve fit the 
actual simulated wetted patterns are shown intable 3 . The points a- 
fone the boundary of the wetted pattern are represented in rectangular 
coordinates, with the same x coordinate for the purpose of statistical 
comparison. The goodness of fit of these points derived from the simu- 
lated model and those obtained from the equations of the ellipse are 
represented by the correlation coefficient that is in the order of 0.989. 
Under no wind condition, the correlation coefficients of circular wetted 
patterns are best fitted as they are equal tol. The possibility of 
using a more sophisticated curve fitting technique to obtain the best 
fit conic equation for distorted wetted pattern is beyond the scope 
of this study. Similar representation of circular wetted patterns were 
used by other researchers (11). 

The distributions of application rates within the boundaries of a 
wetted pattern were assumed to approximate a triangular and semi-el- 
liptical profile as suggested by Bittenger and Longenbaugh (4). Under 
Windy conditions the point of maximum concentration is distorted for 
a distance (c) from the location sprinkler. For a wind velocity of 
20 mph, the distortions for a range of particle sizes 2,3,4, and 5 mm 
are shown in figure 9. The relative displacement of the distortion 
of the centre of maximum concentration for a given wind velocity is 


small (in the vicinity of 2 ft), relative to the throw. The previous 
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distribution profiles can thus be assumed to hold in the distorted case 
with the point of maximum concentration located at the displacement (c) 


of the largest particle (5mm) under consideration. 


6.3 Validity of Interpolation Technique 

Interpolation from the basic wetted patterns were necessarily in- 
corporated in the model for reasons discussed in section 6.1. ‘The 
method of interpolation for unequal space data is noted in section 3.2, 
Chapter 3. Plots of the input variables (d), (IV) and (W) against the 
output parameters (g); (b), (c), (e) and (x), as shown in figures (10), 
(11), (12), (13), (14) and (15), indicated a linear dependence rela- 
tionship exists. However, the unequal slopes of the plotted lines 
suggested the use of a two dimensional interpolation technique (the 
natural bicubic spline approximation). The accuracy of the two di- 
mension interpolation technique was substantiated when correlation co- 
efficients of 0.98 were obtained when interpolated values were compared 
to those simulated (see table 4 ). More than two dimension interpola- 
tion techniques may be used to increase the accuracy of the desired 
output parameters for the purpose of reducing the computer simulation 
time. Accessibility to these programs was not possible during the 


course of this study. 


6.4 Centre Pivot Performance Analysis 

In evaluating the performance of the simulation model eleven com- 
parisons of simulated and field data were available for the two C.P. 
systems tested. The effect of the wind was incorporated in the model as a 
vector for a fifteen minute sequential interval. The triangular and el- 


liptical distribution profiles were used to represent the actual case. 
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Figure 10. Relationship between the semiaxis (a) and the 
input variables. 
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Figure 13. 
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Relationship between the offset from sprinkler (c) 
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Relationship between the semiaxis (b) and the 
input variables. 
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input variables. 


2 S fal: 

20 
A, 
= 
ae 
=a 10 
iS Iv=90 ft/s 
“a W.= 0 
E d 
v 2,3,4,5=drop g (mm) 
@) 
S 
4 
= 


0 20 40 60 80 100 


Semiaxis (a), Feet 


Figure 15. Relationship between the semiaxis (a) and the 
input variables. 


a» eee 


ta’ (4 7 


a of rele Qatar otf eke 


63 


(4) 
JUETOTFJOOD 
UOTIeTSAIOD 


(O)e 
peqnduop 
STxy JoleW 


peyepTodizs3ut 
STxy JoLey 


(AT) 
AATOOTEA 
qetTdoiagd [TetITUT 


(4) 
FUSTOTF FOOD 
UOTILTSEAIOD 


(2)q 
peqndwop 


STXY AOUTW 


(S)q 
poezepTod1zsquy 
STXY JOUTW 


AV TOOTOA 
eTdorg TeTITUT 


“SNNELEVd CELLAM UTTININdS ONINTAHG SUALAWVYVd GALWIOGNALNI CNW GaLVINWIS HHL lO NOSTYVdWOO 


‘Ul «SZTS STOTUWeg 


x 
(-_0T * 33) 


“7 TIUML 


64 


Evaluation of the model involved the comparison of the application 


depths and the coefficient of uniformity. 


6.4.1 Comparison of Field and Simulated Application Depths 

The adequacy of the model for predicting the depth of application was 
determined by comparing with experimental measurements from the two 
operational C.P. systems. Field application depths were measured in 
receptacles spaced 20 ft apart and extending radially outward from the 
center-point (see figure 7a and 7b, section 5.2). The individual 
sprinkler wetted pattern, the angular velocity of the lateral and the 
Wind characteristics were measured and input in the model to allow a direct 
comparison of field and theoretical simulated depths. Attempts to use 
a statistical test for comparison to a standard or control set of field 
measured application depths were not possible because of the chaotic na- 
ture of the influencing environmental factors. The disparity between 
the measured and simulated depths were situated in terms of the absolute 
deviation (AAD1) and percentage absolute deviation AAD2, computed by 


the following relations: 


lO) > Tay! 


AAD]. = 3 (ce) 
[yee eal aul ( 37 ) 
0 
VAD (1) 
N 
where, 


0. - individual measured application depths (ins) 
T. - individual simulated application depths (ins) 


N - number of sets of measurements used in comparison (integer) 
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The theoretical application depths were determined for both tri- 
angular and elliptical distribution patterns and plotted against the 
field measurements that were determined for 20-ft radial intervals. 

Only small differences were observed between computed application depths 
for triangular and elliptical sprinkler patterns (see figures (IV(7) to 
(IV(77) in Appendix Iv). The close spacings of the adjacent sprinklers 
along the lateral probably concealed the effect of the two assumed dis- 
tribution profiles. A similar observation was made by Heerman and 

Hein (11). The triangular distribution profiles were selected for fur- 
ther comparative analysis because they gave a slightly better estimate 
of the actual field application depths. In general, the average abso- 
lute deviations of the theoretical application depths for triangular 
distribution profiles were Ted more frequently than those of the ellip- 
tical profiles (see table (5), (6), (7) and (8). 

Further evaluation of the model includes the comparison of measured 
application depths and those theoretically derived using the triangular 
distribution profile for wind and no wind conditions. Plots of the field 
and simulated depths for all the tests carried out on the two C.P. sys- 
tems are shown in figures (16) to (26). From figures (16) to (20) the 
following observations may be summarized as follows: 

(i) There is generally good agreement between field and theoret- 
ical application depths along the lateral. This is shown by 
the approximately similar trend in the plots. 

(ii) A closer agreement is obtained when the wind magnitude and 

direction is considered in the simulation model. Im this 
case the model attempts to approximate the effect of the wind 
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cases, however, the effect of the wind is minimal (see fig- 

ure 18, 21,26). Either the effect of the wind is insignif- 
icant or an overall "averaging out" effect occurs in the a- 

bove cases. 

(i111) In most cases a large application depth is noted near the 
center of rotation. The application rate near the pivot 
point is low but the time of application is long. As a de- 
Sign attempt to achieve a uniform depth the time of applica- 
tion on a point decreases with a corresponding increase in 
in application rate outwardly from the centre pivot. 

(iv) Fluctuations of the measured depths along most of the lat- 
eral is evident. These fluctuations may be attributed to 
the localized chaotic effect of the wind on the superimposed 
adjacent patterns. The overlapping sprinkler effect may be 
magnified or greatly reduced because of a rapid change in 
wind direction. 

(v) A lag between the field and theoretical depths is noticeable 
near the centre pivot. Inaccurate time synchronization 
in the input of the simulation model is most likely to cause 
these shifts. 

The simulation model is best evaluated in terms of the percentage 
absolute deviation (AAD2) which was computed earlier from equation 37, 
section 6.4.1. On the whole, theoretical depths deviated approximately 
14 percent from the actual field measurements. The spread of average 
deviations for the "Circle Master" C.P. system varied from 11 to 18 per- 
cent as compared to the "Valley 1160" system of 14.3 to 30.8 percent 


(see table 5,6,7 and 8). In general a better agreement was obtained 
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for the former C.P. system as compared to the latter system. A mean devi- 
ation of 14.7 percent and a standard deviation of 1.96 percent in con- 
trast to 19.4 percent and 4.7 percent, respectively, was observed. This 
disparity may be attributed to the individual system performance character- 
istics and statistical differences probably because eight fewer tets were 
carried out on the "Valley 1160" C.P. system. When the effect of the wind 
was not included in the model the mean absolute deviations were increased 
to 17.5 percent for the "Circle Master" C.P. system compared to 19.4 per- 
cent of the "Valley 1160" system. The influence of the wind is not prev- 
alent in the latter C.P. system because of the wide fluctuations of the 
measured depths along the lateral, (see figures 25 and 26). In terms of 
mean absolute deviation, no valid assessment of the theoretical values 


were possible in this case. 


6.4.2 Factors Affecting Model Performance 

The differences in measured application depths and those theoretically 
derived by the simulation model could be attributed to inaccuracies in 
field measurements of significant parameters and variables, changes in sys- 
tem characteristics and in the assumptions used in oe mathematical model. 

Inaccuracies in field measurements of application depths may be the 
result of a number of factors that often occurred in the field. Detri- 
mental factors such as the occurrence of tipped receptacles under wind 
conditions, the influence of surrounding vegetation, malfunction of sprink- 
lers, splashing effect of water streams emitted from leaks, and the prop- 
elling mechanism (Valley 1160) often were not compensated in this study. 
The variation of system discharge due to change in pressure as well as the 


intermittent change in sectional angular velocity of the lateral, are 
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probable causes of the large fluctuation in the measured application 
depths. Difficulty in measuring the chaotic nature of wind magnitude 
and direction of the localized area is another probable reason for 

these fluctuations in field measured depths. The overlapping effect of 
the sprinklers is either enhanced or reduced depending on the wind con- 
ditions at a particular given tine. 

Errors in the theoretically determined depths are the rusult of in- 
accurate input variables and the assumptions made in the model. Proper 
synchronization of the starting and ending time of the field experiment 
is necessary to determine the appropriate range of monitored wind input 
data for the simulation model. The length of experimental time often 
was not easily estimated for rain gauge monitors. In the simulation 
models, the effect of the wind is considered in sequential times of 15 
minutes and the chaotic nature of the actual wind conditions were consid- 
ered aS a one dimensional vector. The error incurred was necessary be- 
cause of the long computer simulation time of the test run as well as 
the limitation of wind monitoring equipment. Therefore, the influence 
of the localized chaotic windy condition on the application depths was 
not accurately modelled in the program. The range of particle size 
estimated using the method discussed in Section 5.4, Chapter 5, is only 
an approximate estimate. If the distributions of the particle sizes 
for each nozzle were specified the desired wetted patterns would then 
be very accurately determined for the flight trajectories. 

The distorting effect of the wind on the assumed triangular and 
elliptical distribution profiles were not considered in the simulation, 


since the particle size distribution pattern was not known. However, 
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an attempt was made to determine the point of maximum concentration 

when the patterns were distorted (see section 6.2). This point of max- 
imum application depth was idealized asthe origin of the distorted wetted 
pattern. Accurate determination of the point of maximum concentration 
made under the chaotic influence of wind was not possible. Error in 
assuming a triangular or elliptical distribution profile is relatively 
small (11), however, the distorting influence of the chaotic wind con- 
ditions could greatly emphasize or reduce the overlapping effect of the 
profiles of adjacent sprinklers. 

The effect of a change of wind direction on the application depth 
could be theoretically studied by using a constant wind direction rela- 
tive to the pivot throughout the analysis. For wind velocity of 0, 20, 
40 mph and a change of direction of 0, 90, and 180 degrees, relative 
to the pivot, the application depths for both C.P. systems are shown 


in Appendix V, figures Vay Cony The effect of constant wind direc- 


(G)s 
tion was minimal throughout for most of the lateral. Disturbance at 


the beginning and at the end was most obvious when wind direction is at 
0 degree relative to the pivot. The plots for windy conditions (20, 


40 mph) as listed in figures Via) and Via)? showed a shift of the curves 


towards the centre of the pivot. For a wind direction 180 degrees rel- 


ative to the pivot the opposite effect was observed (see figures V3) 


and V . An intermediate response of these two extremes was observed 


(6) 


for a wind direction of 90 degrees in the figures V and Ves. ihe 


(1) (5) 
effect of change in wind velocity was insignificant in the application 
depths for constant wind direction. It may be noted that the effect 


of chaotically changing wind direction has a greater overall effect on 


the distribution of the application depths along the lateral than 


Ci ai 
1a puis 


cada bed base tt $8 orld 6B» 
ts ge a ‘y we 
Poet ot et ee) rence Ser ten wth iy 
ut Sore: in 


@ ie 
Wisviis lai ae ot Lies raped. 


“nce *ciw oloeaa “all ng ponesC ght 


a to. PaaTts cart gee Pope iui autor ach 


efey qnicbet th Sale 24 » priee a tai cae han 

Oe fo ealoetler Biew 1 acer a ag ees ne ty ests os 

wtktntes .enseoek mie Ue 40 at reat sae? ate teacris 6s bt a OF 

pworie Sw ero. 2 big ek ache eed feces, ene strove mit oF 

eet Ey iwi JT HK Ww wae ergz 4 ‘ ai Oy Ty VY. ie td ; oi eS 
‘ ie 

: eee fab 

& evi ner seo te Meee 2c! cig ee Penn 7 

& TT Pare Th ; rh, wa, ean «4 A ors) owt 77 aS acini - y. 

arottihrye vont Sel seh a baby weld Oo oebeatert - snes 

mH ow Tar Ca v Sage) ri Soe ni iaiahaan Tee 


if pepe a 
Peay : Sx ang Phy) 2 Yoda Se ae Bi oe ait) i sires iain hd 4 


ay: =) Lengel: ate spy ara seal orkg was oe a koe 
i ; : a 7 See 


hevabrt: tes. eect ce ser’ bo eX eed 2) 


MEE ce ng WR ihe 4 awk he orl oe seat | 


Pach a aaitt ok Hn sansiclaale a vate 
ais “OPEN stage Pedtion a, eee Hoes ls up te re " * 


me it 
ibe , Ai Ora ae ce 


ehh esate are it ts baa pceniert, 


85 


the magnitude of wind. 


6.4.3 Coefficient of Uniformity 

Christiansen's coefficient of uniformity was used to evaluate the 
performance of the C.P. systems. A similar weighted coefficient of 
uniformity suggested by Heerman and Hein (11) is also included in this 
analysis to give more weight to the application depths further from 
the pivot where the area of influence is increasingly larger. 

Comparison of the field and model simulated coefficients for both 
weighted and non-weighted values are shown in tables 5,6,7 and 8. The 
uniformity coefficients in this table were calculated by excluding the 
area at the outer radius where the large nozzle (end gun) sprinklers 
are located. The deviation of the mean of the sprinkler system is af- 
fected by a large area at the end guns and reduces the coefficient ac- 
cordingly. The coefficients of uniformity determined for both triang- 
ular and elliptical distribution did not differ greatly (see table 9). 
As discussed earlier (see section 6.4.1) the overlapping effect of 
closely spaced adjacent sprinklers for the two assumed distributions 
gave only small differences. In general, however, the theoretical sim- 
ulated coefficients were higher than those measured. Thisis evident 
because of the wide fluctuations in field measured depths and the in- 
itial larger response close to the pivot. Closer agreements were ob- 
tained for the measured weighted coefficients and those simulated in the 
model as a standard deviation of 1.26 compared to 1.56 was noted in the 
above case. The accuracy of the simulation model cannot ke definitley 
assessed irrespective of the reasonable agreements of the theoretical 
computed and measured coefficients. The coefficient of the uniformity is 


only a performance parameter of the application depths along the lat- 
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eral. The variation of application depths at any given point beneath 
the lateral cannot be fully accounted for without imposing an averaging 
effect. The effect of wind conditions on the coefficient of uniformity 
is theoretically predicted in table 10 and 1l. In this analysis the wind 
direction was held constant is in a direction relative to the centre 
pivot. At 0 degree direction where the line of measurement is directly 
Opposite to the wind a noticeable reduction in uniformity coefficient 
occurs. The extreme distortion of the overlapping distribution patterns 
in this direction accounts for the lower coefficients. From tables 9 


and 10 the coefficients obtained under the range of wind magnitude and 


direction are within the acceptable 80 to 90 percent design requirements. 


6.4.4 C.P. Field Performance Analysis 


A summary of the field measurements that were conducted (as dis- 
cussed in Section 5.1 to 5.4, Chapter 5) are shown in table 12 and 13. The 
average application and coefficients of uniformity were computed for the 
various test runs in each C.P. system. In each C.P. system the average 
application depth and coefficient of uniformity showed only small vari- 
ation ( Std. Dev. of approximately 1.1) for each selected speed of 
lateral rotation. Although a wide range of wind velocities were meas- 
ured , the average wind directions that were recorded did not indicate 
the actual. fluctuations recorded during the experimental test period. 

The possibility of determining the performance of the two C.P. 
systems in relation to the varying environmental parameters from the 
results shown in table 12 and 13 was’ restricted. The available meas- 
urements of average wind directions and velocities in tables 12 and 13 
cannot be related to the determined application depths and coefficients 


of uniformity because of long experimental runs of approximately 14 
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TABLE 10. PREDICTED COEFFICIENTS OF UNIFORMITY FOR THE "CIRCLE MASTER" 


SYSTEM. 
Triangular Elliptical 
Distribution Distribution 
Wind Wind Pattern Pattem 
Velocity Direction 
M.D iN. 3 g % % 
0 95.82 95.00 94.86 93.33 
0 87.30 9205 86.91 91.76 
50 90 95.56 91.72 94.85 90.89 
180 96.02 89.88 95/2 89.49 
0 Apes S528 74.14 85.90 
40 90 89.35 81.25 89.80 81.40 
180 93.46 82.07 93.37 81.28 


TABLE 11. PREDICTED COEFFICIENTS OF UNIFORMITY FOR THE "VALLEY 1160" 


S¥STEM. 
Triangular Elliptical 
Distribution Distribution 
Pattern Pattem 
Wind Wind 
Velocity Direction Cus Cu. Gus ese 
(Wt.) (N.Wt.) (Wt. ) (N.Wt.) 
m.p.h. % % g % 
0 95.96 OSS Say sul 93.59 
90.81 
94.29 
04250 
0 78293 88.58 79.56 88.60 
40 90 95.30 86.82 93.56 86.07 


180 9500 Soo 94.00 Go.02 
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hours involved in each test. The measured field data, however, provided 


a guide in the development of the C.P. simulation model. 
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CHAPTER VII 
SUMMARY AND CONCLUSIONS 

The main objective of this study was to evaluate the performance 
of two Centre Pivot irrigation systems operating under local climatic 
conditions. A mathematical model was developed to simulate the dynamic 
operation of a C.P. system. The simulation model could then be used as 
a tool to determine the system performance under actual windy conditions. 

In the first part of the model the dynamics of water droplets in 
flight were formulated as second order differential equations to deter- 
mine the basic wetted patterns. The required input information such as 
che range of droplet particle size and initial particle velocities were 
obtained from experimental analyses and the manufacturer's sprinkler 
specification, respectively. Basic information on wind direction and 
Magnitude in discrete time intervals of 15 minutes were considered in the 
analysis to study the distortion of the wetted patterns and the overlap- 
ping effect of the adjacent sprinklers. The wetted pattern of sprink- 
lers for no wind condition is described by an equation of a circle. Un- 
der windy conditions, however, the wetted pattern is approximated by 
curve fitting using two overlapping ellipses with common minor axis. 
The point of maximum concentration in this case is distorted a distance 
(c) from the origin (i.e., the sprinkler's location) . 

Triangular and elliptical cross-sectional distribution profiles 
were assumed to represent the extremes for individual sprinkler patterns. 
Interpolation techniques were used to determine the desired wetted pat- 
terns from those derived from closely related initial conditions. 

This was necessary because of the extremely long computer simulation 


time required to compute the actual wetted patterns from the droplet 
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trajectories. 

The second part of the model involved the simulation of the oper- 
ating C.P. system in relation to the rows of measuring receptacles used 
in field measurements. Actual simulation of the moving wetted patterns 
for the rotating sprinklers were developed for % degree incremental ang- 
ular change of the lateral. Actual wind information monitored during 
the testing period were modelled in the program to assess the perform- 
ance of the simulation model in terms of application depths and the co- 
efficient of uniformity. 

From the development of the performance from available experimental 
data the following conclusions were drawn: 
iu Wetted patterns of a rotating sprinkler on a moving C.P. system 

can be defined from the equation of motion of sprinkler droplets 
emitted. 

Pay The distorted wetted patterns (under windy conditions): can be ap- 
proximated using the equations of two ellipses having a common 
minor axis. The point of maximum concentration of particles in 
this case is a function of the magnitude and direction of the 
wind. 

ci The range of particle sizes that was obtained experimentally for 
use in the simulation model was an approximation. For greater 
precision in the application of the simulation model more accur- 
ate techniques are suggested. The differences between the ac- 
tual wetted patterns and those simulated are most likely greatly 
affected by the discrepancies in the selection of appropriate 
particle sizes. Range of particle size is the significant cri- 


teria in determining boundaries of the wetted patterns. 
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The interpolation techniques used such as bi-cubic spline were 
valid in determining the desired wetted patterns for the purpose 

of reducing the extremely long computing time. Correlation co- 
efficients of 0.98 were obtained for interpolated and simulated 
parameters describing the wetted patterns. 

Comparisons of experimental and simulated application depths for 
most test runs gave good agreement. A deviation of approximately 14 
percent of the simulated , compared to that of measured, was ob- 
served. Wider fluctuations occurred in the field. This is attrib- 
uted to the chaotic nature of the wind and discrepancies in field 
measurement techniques. However, a reasonably similar trend was 
obtained in the plot of application depths simulated by the com- 
puter. 

Triangular and elliptical distribution profiles that were assumed 
in the simulation model did not differ under most conditions test- 
ed. In general, however, the triangular distribution profile 

gave slightly better agreement with the actual distribution profile. 
Time synchronization of the experimental process and accuracy in 
recording measurements of influencing environmental parameters 
(such as wind velocity and direction) is required to improve the 
accuracy of the model so that the wide fluctuations in the field 
are accounted for. The simulation model can then be accurately 
validated. 

The two Christiansen coefficients of uniformity used in this 
study gave only a satisfactory measure of uniformity. A better 
agreement is obtained when: the uniformity is evaluated using 


the weighted coefficient. In general, although wide fluctuation 


i 


; ae . paris at 


| 
beet 2: hn) Beis 
’ 5 od 4 tf 
per mT 
22 abtelengrne th tm, ach ole SP pitti ais 
plimke “Sous an : deceleebrioes scones 
Ae 
. | wots Joie elt | 
: ein A Grace ST: z : 
a z dy, i wen pe- i Ly } | fe i 
f ie ) 4 n 
) ae 
_ 
ie SiS EY ; eben, es ; 
\ . hy ‘9 sas tical 7 
[ Fi ty et Oke ash enue a 
an 
ay stim Pople oF 
ae. ee “ut eet P ae | 3 ae “a 
Ve a 
sith’ SA ca RRS Es 23% Ftiemops oan 
| lw, to eet Howe, Sih ste 
] il 7 
rasiine de <laeeee emioms sale AS _ ae 
: >> Peper Pies oy rerio ‘ban : au ia bs Petey 
as Shih Ts ote Alien 2a 
y - 


1a 


in application depths occurred under the experimental conditions , 
the coefficients of uniformity derived fell within the higher lim- 
its of acceptance. The overlapping effect of the closely spaced 
adjacent sprinklers is most probably resposible for such values 
of coefficient of uniformity. 

As a result of a number of field experimental errors in the 
measurements of application depths, environmental parameters and 
long experimental test times, no valid relationships can be made 
between the coefficient of uniformity and the variation in the 
environmental factors that occurred during the experiment. 

The simulation model provided a method for studying the effect of 
Wind magnitude and direction which was a continuous phenomenon 
throughout the experimental period. With proper validation of 
the simulation model, further improvement in design criteria for 
C.P. sprinkler irrigation systems can be made. For example, the 
optimum spacing of sprinklers, sprinkler type, and size of the 
nozzles for a given set of known system characteristics, can be 


theoretically determined. 
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CHAPTER VIII 


SUGGESTIONS FOR FURTHER STUDY 


Further research oriented toward improving the accuracy and the 


efficiency of the simulated model should consider the following areas 


of work: 
i By More accurate determination of the sprinkler droplet sizes. 
Ze Construction of the sprinkler distribution patterns entirely on 


the basis of droplet size distribution within the sprinkler 
wetted patterns. 

Se Implementing losses due to evaporation and wind drift into the 
equations of motion. 

4. Elimination of the interpolation procedures used in the model. 

5 Validation of the simulated model based firstly on experimental 
research in the laboratory under controlled climatic conditions 


and secondly under actual field conditions. 
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APPENDIX II. FLOW CHART OF THE OPERATING SIMULATION MODEL. 


Read wetted 
Dbaecerns, 
system 
baka. 


: S 
Doglinear: 
Interpolation 


Continue 


Input variables defining 
basic wetted pattems and 
system specifications. 


Input position of the lateral, 
direction of wind, wind velocity, 
number of 4° increments of 
lateral rotation. 


Do for each sprinkler. 


Do two dimensional 
interpolation. 


If wind velocity greater 
than zero, do linear 
interpolation. 


Set counter to zero, 

Count nosiof #5 incr. of at. rota. 
Generate position of lateral. 
Generate constant wind direction 
for given no. of k°®° increments. 
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APPENDIX II. CONTINUED 


Compute application rates and appli- 


cation depth for each point of the 
first row.* 


Compute application rates and appli- 
cation depths for each point of the 
second row.* 


Compute Compute application rates and appli- 
TR3) His cation depths for each point of the 
AWAD Ses abies: third: rowe* 


>) [continue | 
Write application depths for each 
of the three rows. 


*See detail flow chart on the following page. 
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Do for each 
pointe (P). 


Do for each 
sprinkler (Q). 


Consider: 
Circular patterns (W=0) 
or Ellip. patterns (W=0). 


Compute distance from 
point of max. appl. rate 
(Q or A) to point under 
consideration (P) and 
toj poine, on, badry..") (Z).. 


Seakeh Li pera) is 
Within the wetted pattern. 


Compute appl. rate and 
partial appl. depth for 
Cr lang sono celllpce.c1str. 


Search if any other 
wetted pattern encloses 
DOIN. (Pye 


Search if any other 
point (P) receives 
water application 


Compute total appl. 
depth for each point (P). 
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APPENDIX IV 


This appendix contains figures IV to figures IV 


(1) (11) a 
showing the comparison of field and simulated application depths 


for the "Circle Master" and "Valley 1160" C.P. systems. 
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APPENDIX V 


This appendix contains figures Vig 22 Vig)! showing the 


distribution of application depths for the "Circle Master" and 


"Valley 1160" C.P. systems. 
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